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ABSTRACT 

We present the results of the first theoretical non-radial non-adiabatic pulsational 
study of M dwarf stellar models with masses in the range 0.1 to 0.5 M0. We find the 
fundamental radial mode to be unstable due to an e mechanism, caused by deuterium 
(D-) burning for the young 0.1 and 0.2 M© models, by non-equilibrium Hc'^ burning 
for the 0.2 and 0.25 M0 models of 10''' Myr, and by a flux blocking mechanism for the 
partially convective 0.4 and 0.5 M0 models once they reach the age of 500 Myr. The 
periods of the overstable modes excited by the D-burning are in the range 4.2 to 5.2 h 
for the 0.1 M0 models and is of order 8.4 h for the 0.2 M0 models. The periods of the 
modes excited by He'^ burning and flux blocking are in the range 23 to 40 min. The 
more massive and oldest models are more promising for the observational detection of 
pulsations, as their ratio of instability e-folding time to age is more favourable. 
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1 INTRODUCTION 

M dwarf stars are cold, low-mass stars on the main sequence 
that have recently gained attention a s being the best can - 
didates to harbour Earth-like planets jPelfosse et al.|[2006l ). 
Although they make up about 75% of the stars in our Galaxy 
and the solar neighbourhood, their main physical parame- 
ters, such as their radii, are poorly understood, and hard to 
determine observationally due to their low luminosities. 

Evidence has been accumulating for discrepancies be- 
tween theoretically predicted radii and effective tempera- 
tures of low mass stars (LMS), and those derived from obser- 
vational data. Observational studies of LMS in eclipsing bi- 
naries and single LMS have shown that their radii are about 
4 to 20% larger, and their temperatures as mu ch as 5% lower 
than predicted by theoretical models (see e.g. iMorales et al] 
I2OO9I and references therein) . This discrepancy is usually at- 
tributed to effects a ssociated with chrom ospheric activity or 
metallicity (see e.g. lL6pez-Moraie3l2007l ). 

Another reported problem is age discrepancy: although 
mean ages derived from lithium isochrones for LMS in young 
groups agree with those derived from Hertzsprung-Russell 
diagram (HRD) isochrones, the two methods often give dis- 
cordant ages when a pplied to the indivi dual stars in the 
young groups (see e.g. lYee fc Jensenll20ld -YJIO). 

Magnetic inhibition of convection provides a possible so- 
lution to both problems. This leads to models having larger 
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radii and lower temperatures (iMuUan fc MacDonaldl 2001 ). 
and c an remove the age discordance ( MacDonald fc MullanI 
l2010l ). at least for the stars in the YJIO sample. A different 
solution to the radius problem is given by IChabrier et al.l 
(2007) who reduce the mixing length parameter (a) which 
also gives models with larger radii and lower tempera- 
tures. The two approaches have a different signature on 
the stellar models: while reducing a affects mainly the 
outer super-adiabatic regions of a surface convection zone, 
magnetic inhibition of convection affects the whole of the 
convection zone. This produces different profiles for the 
sound speed (cs) in the interior of the star, and raises the 
possibility of differentiating between the two approaches 
through asteroseismic techni ques. A preliminary calculation 
iMacDonald fc MullanI bOlOD using a magnetic and a non- 
magnetic stellar model show that the so called large and 
small frequency separations (related to Cs and its derivative, 
respectively) given by the asymptotic theory iTas soul 198^) 
could distinguish between the two theoretical formulations, 
reaching differences in the two quantities up to 2 and 30%, 
respectively. 



Here we present the first non-radial, non-adiabatic pul- 
sational study of M dwarf stellar models. The destabilization 
of the fundamental radial mode is achieved through an ac- 
tive e mechanism caused either by deuterium (D-) or He'^ 
nuclear burning. For the more massive and older models 
considered, the blocking of the radiative flux at the base of 
the outer convective envelope explains the destabilization. 
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2 PREVIOUS WORK ON LOW-MASS STAR 
INSTABILITY 

The excitation mechanisms known to operate in solar hke 
and LMS are the stochastic and, also for the latter, the 
e mechanism. Stochastic oscillations are driven by turbu- 
lent convection and have been detected in solar-like, j3 
Cephei and red giant stars, which have large outer con- 
vection zones. The t mechanism is produced by the high 
sensitivity of certain nuclear reaction rates to temperature, 
and always has a destabilizing effect. However, no obser- 
vational evidence has been reported so far of stars pulsat- 
ing due to this mecha nism, despite theoretical predictions. 
IPalla fc Bara^ (|2005l ) found theoretical excitation of the 
fundamental radial mode in very low mass stars and brown 
dwarfs (0.02 M/M0 ^ 0.1) due to D-b urning. A non - radia l 
instability study of similar models by iMova et al.l (|201ll ) 
found the excitation of non-radial modes due to the same 
mechanism. Observational varia bility of bo th type of ob- 
jects has been reported (see e.g. ICodvl [20091 and references 
therein), although it has not yet been unequivocally at- 
tributed to pulsations. Photometric variabil i ty of M dwarfs 
toavenport et aLll2011 Hartman et al.|[201ll . iBecker et al.l 
I2OIII . iKowalski et aLll2009( ) has been usually attributed 
to flares and chromospheric or atmospheric activity. Re- 
cently, the results of an on-going time-series photometric 
campaign specifica lly searching for pu lsating M dwarfs has 
been published bv lBaran et al. (|201lh with no positive de- 
tections. 

Previous theoretical stu dies abou t M dwarfs instability 
were done by several authors. iGabriell (|l964l . [l"967l ) proposed 
for the first time the instability of completely convective 
stars on the main sequence, based on the evaluation of a 
vibrational stability coefficient for 0.27 and 0.16 M0 models. 
The fundamental radial mode and its first two harmonics 
were excited when the effects of convection were considered. 

iTomal (|l972h predicted the instability of the fundamen- 
tal radial mode for models of 0.2, . 6, 1.0 and 2.0 M0 due 
to D-burning; while iBourv fc Noels! (|l973l ) found instability 
of 0.5 M0 models, with effective temperatures > 4000 K, 
due to H-burning reactions. 

lOpoien fc GrossmanI (| 19741 ) did a pulsation study of 
models from 0.085 to 0.5 M0 on the H main sequence and of 
0.012 to 0.5 M0 on the D-main sequence. They investigated 
only the fundamental radial mode, which they found to be 
excited due to an e mechanism of the corresponding domi- 
nant nuclear species. They found pulsation periods between 
10 and 41 min for the H-burning models and between 1 to 
20 h for the D-burning ones. Due to the long e folding times, 
they concluded that the mode ls were marginal can d idate s 
to show detectable pulsations. ICabriel fc GrossmanI (119771 ) 
repeated this work with improved physics in the models ob- 
taining pulsation periods in agreement with those previously 
calculated and e-folding times lowered up to a factor of 4. 
The authors pointed out that in order to assure the insta- 
bility of the more massive models, the full non-adiabatic 
equations should be solved. 

The work we present here solves for the first time for 
M dwarf models the complete non-adiabatic pulsation equa- 
tions for radial and non-radial modes. All the previous works 
in the literature are adiabatic oscillation studies, some of 



Table 1. Main pliysical parameters of the excited models. The 
period (Pq) and c-folding time {rfaid) of the fundamental mode 
(£=0, fc=0) are also given. The last column gives e^*/!*^^! for the 
models excited by D-burning. See text for details. 
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them only radial, which based their results in the evaluation 
of a vibrational stability coefficient. 



3 EVOLUTIONARY MODELS 

We h ave used MacDonald's fully implicit stellar evolution 
code ([MacDonald fc Mullan|[201Gl and references therein) to 
construct evolutionary sequences of model stars of compo- 
sition X = 0.6958, Y = 0.2855, Z = 0.0187, and masses 0.1, 
0.2 0.25, 0.3, 0.4, and 0.5 M0. Each sequence is begun at a 
point on the Hayashi contraction phase before the start of D- 
burning. Bohm-Vitense convection is used with a = 1.699. 
The reaction network follows the detailed evolution of D- 
and He"^, and includes all reactions in the ppl chain as well 
as the dominant reactions of the ppH and ppIII chains. 



4 PULSATION ANALYSIS 

We have used the non-r adial non-adiabatic osci l lations code 
GraC o (Granada Code. iMova fc Garridoll2008l . [Mova et all 
l2004h to search for pulsations in our M dwarf evolution- 
ary models with ages from prior to the onset of D-burning 
and up to the age of the Universe. We calculated radial 
and non-radial modes from £=0 to 1=3 in the frequency 
range 20 ^Hz up to the acoustical cut-off frequency for 
an isothermal atmosphere of each model, given by 
I'c = Cs/4:TtHp, where Cs is the sound speed in the sur- 
face of_th^_star_aiid_Hpth^ pressure scale height 
(see iRodriguez-Lopez et al.l I2OIOI for a description of the 
method). 

Our results indicate that only the fundamental radial 
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Figure 1. Differential (dW/dr) and total work (W) done over one oscillation cycle for the 0.1 M0 model at the peak of maximum 
D-burning (age 1 Myr). Also shown are the radial displacement {^r/r), the opacity (logK) and the nuclear energy generation (ejv). All 
quantities are normalized to unity except ejy. Right; Zoom of the surface of the model (r/R = [0.99, 1.0]) 



mode is excited due to an e mechanism caused by the nu- 
clear D-burning (H^(p,7)He'') operating in the young 0.1 
and 0.2 M0 models; the non-equilibrium He''- He"' reactions 
(ppl and ppll chains) operating in the 0.2 and 0.25 M0 
models with age 10* Myr; and a flux blocking mechanism in 
the partially convective 0.4 and 0.5 M0 models older than 
500 Myr. For the radial pulsations excited by D-burning, the 
periods, which are related to the dynamical time scale of the 
object, Tdyn ~ (Gp)^^^^ , are found to be 4.2 to 5.2 h for the 
0.1 M0 and 8.4 h for the 0.2 M0 models. The periods of 
the modes excited due to He'' burning and the flux blocking 
mechanism are in the range 23 to 40 min. 

Table [T] gives the age and main physical parameters 
of the excited models and the period and e-folding time 
of the fundamental mode. Given that the amplitude of 
the mode cannot be calculated with a linear oscil- 
lations code, we refer to the e-folding time, defined 
as Tfoid = 1/|(T/|, where aj is the imaginary part of 
the eigenfrequency, as an estimation of the observ- 
ability of the mode. As the time-dependence of the 
mode amplitude is given by Ae"'^ , where A is the 
unknown initial amplitude of the mode, the e-folding 
time, Tfoid, gives the time needed to increase the ini- 
tial unknown amplitude a factor of e. Therefore, e- 
folding times shorter than the age of the model, or 
than the time scale of the nuclear reaction, would 
allow enough time for the amplitude of the mode to 
be developed, even if the initial amplitude were low. 
If it were large, then even short e-folding time scales 
could allow for observable amplitudes. 

The D-burning lifetime for the 0.1 M0 models is 
1.4 Myr, defining the D-burning phase as at least 50% of 
the star's luminosity provided by D-burning. With Tfoid of 
the order, or longer, than the age of the model, the 0.1 M0 
models are not favoured for the observational detection of 
pulsations, unless their initial amplitudes are large. The last 
column of Table[T]gives an estimation of the possibility of the 
mode being observed th rough the quantity e'^'^''^^' (|Tomal 
Il972l . iMova et al.ll201ll ). This gives an estimation of how 
much the mode grows in amplitude during the time that the 
model spends in the D-burning phase, given by At. With a 



maximum value of 3.3 for the 1 Myr model, at which the 
D-burning reaches its maximum, the probability of being 
observed is not very much favoured. 

Figure [T] (left) shows the work integral (W) and its 
derivative (dW/dr) for the fundamental mode of a repre- 
sentative 0.1 M0 model. The epsilon mechanism asso- 
ciated to the D-burning is responsible for the driv- 
ing {dW/dr > 0), as the work integral in the cen- 
tral region of the star is mostly dominated by the 
temperature derivative (et ~ deN / dT\p) of the nu- 
clear energy generation (ejv), which in our models 
has a high value, et — 12 at temperatures T ~ 10^ K. 
The epsilon mechanism is strong enough to over- 
come the damping of the outer regions of the star 
when er is large and the central density is low, al- 
lowing for larger amplitudes (Cr/r) in the inner re- 
gion of the star. As the burning of D- proceeds, er 
decreases from its initial value er — 12, at tempera- 
tures T ~ 10® K, for the excited models to er — 6 for 
the stable ones. Figure [T] (right) zooms in the surface of 
the star to reveal a small damping {dW/dr < 0) and a large 
driving region. However, they do not contribute significantly 
to the total work, as ~90% of the work integral is given by 
the e mechanism in the inner regions. Even if the dW/dr 
profile in the outer layers seems to be produced by 
a K-mechanism, the maximum driving should then 
be produced at maximum opacity, i.e. in Figure [T] 
lognmax at r/R ~ 0.93 or equivalently at logT ~ 4.6, 
which corresponds to the partial ionization zone of 
Hell. However, the maximum driving is produced at 
r/R ~ 0.998 corresponding to logT ~ 3.6, where no 
partial ionization zone is expected. A careful analysis 
of each of the terms in the energy equation shows that is not 
the K mechanism which plays a role in the destabilization, 
but the radiative luminosity, which is very model dependent 
and not trustworthy for these convective or quasi-convective 
models. 

We follow iMova et al.l (|201ll ) to assess that the insta- 
bility of the modes is really produced by the e mechanism: 
checking the invariability of the work integral for different 
mechanical outer boundary conditions. We apply the crite- 
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Figure 2. Normalized work (W, blue dashed lines) and its derivative (dW/dr, blue solid lines) for the fundamental radial mode of 
the 0.5 M0 model of 500 Myr with and without nuclear reactions (ejv=0, red lines). Right: Zoom of the previous plot in the range 
r/R = [0,0.99] to see the details of the driving at the base of the outer convection zone, located at r/R = 0.57 



rion that for instability at least 80% of the work integral 
must be provided by the e mechanism. 

The D-burning lifetime decreases with increasing mass, 
being 0.7 Myr for the 0.2 M0 models. Only the 0.2 M0 
model at the peak of maximum D-burning (0.5 Myr) is ex- 
cited via the associated e mechanism, although with a low 
probability of being detected, as Tfoid is of the order of 
the age of the model and consequently e'^*'^''^^' has a low 
value. The destabilization of the radial fundamental mode 
for 0.2 and 0.25 M0 models of lO'^ Myr is caused by the e 
mechanism powered by the He^ burning produced in the 
pp chains. The period of the modes is 23 and 27 min, re- 
spectively, with e-folding times low enough for the pulsating 
amplitude to be developed, although the observational chal- 
lenge for these elderly objects seems insurmountable. 

The 0.4 and 0.5 M0 models are partially convective: 
they have an outer, and the younger ones also an inner, 
convection zone. The excitation of the fundamental radial 
model for ages older than 500 Myr is due to flux block- 
ing at the b ase of the outer ex ternal convection zone, as 
described by iGuzik et all (|2000l ) for 7 Doradus stars. The 
radiative luminosity falls abruptly at the bottom of 
the convection zone as it can not be transported by 
radiation. For these models, as the local convective 
life time at the base of the convection zone is much 
longer than the pulsation period, the convective flux 
can not adapt instantaneously to transport the ad- 
ditional energy coming from below. Therefore, this 
energy is periodically blocked and converted into 
work driving the oscillations. The large local con- 
vective lifetime compared to the pulsation period 
validates our choice of fr ozen convection approxi- 
mation (|Unno et al.l Il989l ) in the oscillation code, 
where variations of the convective flux are not con- 
sidered during the oscillation. Fig[2](left) shows the to- 
tal work (blue dashed line) and differential work (solid blue 
line) for the fundamental radial mode of the representative 
0.5 M0 model of 500 Myr. The same quantities are shown 
with red lines for the same model with the nuclear reactions 
set to zero. The driving region r/R = [0, 0.57] is due to the e 
mechanism associated to the He'' reactions. When the con- 



tribution of the nuclear reactions is switched off, radiative 
damping is the main contributor to the energy. However, the 
mode is still excited due to the flux blocking mechanism at 
the base of the outer convective region at r/R=0.57. This 
is better seen in FigO (right), where we have left out the 
r/R=[0.99,l] region, in which the effects of the radiative lu- 
minosity on the driving are not trustworthy. 



5 DISCUSSION AND CONCLUSIONS 

We presented the first non-adiabatic non-radial pulsational 
study of M dwarf models. We found the fundamental ra- 
dial mode to be excited by an e mechanism caused by D- 
or He"' burning operating in the fully convective models, 
or by a flux blocking mechanism for the partially convec- 
tive ones. The 0.1 and 0.2 M0 models excited by D-burning 
have growth time scales of the order of the age of the mod- 
els, which makes their observational detection elusive, unless 
their in itial amplitudes are lar ge. This contrasts with the re- 
sults of lPalla fc Baraffd Hooi) for 0.1 M0 models of 0.6 and 
1.1 Myr which have e folding time scales more favourable for 
detection. 

Oscillations in the 0.2 and 0.25 M0 of 10"* Myr, excited 
by the e mechanism caused by the He^-He^ reactions, are 
not envisaged to be detected observationally, due to the age 
and faintness of the models. 

The observational detection of oscillations for the 0.4 
and 0.5 M0 models excited by a flux blocking mechanism 
is more likely, as the amplitude of the modes have enough 
time to develop due to the shorter time scale for the growth 
of the instability. 

We now discuss whether the predicted pulsation 
periods could be misinterpreted in observational 
light curves as rotational periods. IJenkins et al.l 
(|2009h (J09) derive vsini for a sample of 56 M 
dwarfs with spectral types M3-M6, corresponding 
to masses in the range 0.1-0.3 M0, obtaining a dis- 
tribution that peaks at vsini = 3 kms^^. They add 
to their data a sample of about 200 M dwarfs with 
spectral type in the range M0-M9 and vsini calcu- 
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Table 2. Estimation of rotational velocities (2ttR/P) for the low- 
est and highest mass range of the models in our study. 
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lated in the literature. They found a large spread 
in vsini, from about to 50km/s, for the mid spec- 
tral range M4-M6.5 (see figure 1 of J09), although 
for a large fraction of these stars the measured vsini 
are upper limits. In general, there is a trend of in- 
creasing rotational velocity with later spectral type, 
i.e. decreasing temperature and mass. We roughly 
estimate the rotational velocities (2-kR/P) assuming 
that the pulsation period of our star models is the 
rotation period (see Table [2|) . 

We conclude that the pulsation periods of the 
0.1 and 0.2 A/0 models (corresponding to spectral 
types about M5-M9) excited by the epsilon mecha- 
nism (cases a, b and c) can not be confused with 
rotation periods, which have maximum values of 
vsini < 50 kms~^ (figure 9 of J09). For the highest 
mass range of our study (0.4-0.5 A/0 models, cor- 
responding to spectral classes about MG-M2) where 
pulsations are excited by the flux-blocking mecha- 
nism, the pulsation periods would be distinguish- 
able from the rotation periods (vsini < 10 kms~^, 
see again figure 9 of J09). 

Furthermore, a M dwarf star that rotates with 
a period as short as the pulsation periods given in 
Table [l] would likely stand out as a flare star, which 
would make it an unsuitable target to detect pulsa- 
tions. 

The observational confirmation of pulsations in M dwarf 
stars would mean to open the field of asteroseismology to 
characterize this class of stars with much more precision 
than it has been done before, as an independent measure 
of the stars' mean density could be directly made 
from the oscillation periods. This would result in 
a more precise physical characterization of possible 
planets hosted by pulsating M dwarfs. 

Finally, we can not discard that modes of high 
radial order in the asymptotic regime are excited 
by an stochastic mechanism. The study of this pos- 
sibility, with dedicated codes that include non-linear 
interactions in the convection zone, would be desir- 
able as it could lead to resolution of the oversizing 
puzzle for M dwarfs described above. 
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